A repulsive reference potential reproducing the dynamics of a liquid with attractions 
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A well-known result of liquid state theory is that the structure of dense fluids is mainly determined 
by repulsive forces. The WCA potential, which cuts intermolecular potentials at their minima, is 
therefore often used as a reference. However, this gives quite wrong results for the viscous dynamics 
of the Kob-Andersen binary Lennard- Jones liquid [Berthier and Tarjus, Phys. Rev. Lett. 103, 
170601 (2009)]. This paper shows that repulsive inverse-power law potentials provide a useful 
reference for this liquid by reproducing its structure, dynamics, and isochoric heat capacity. 



Condensed matter is held together by attractive inter- 
molecular forces, but is generally very resistant against 
compression. These two simple observations provide im- 
portant microscopic information, and long ago it was 
shown that many properties of liquids and solids derive 
directly from the fact that intermolecular forces consist 
of short-ranged, harsh repulsions and long-ranged, much 
weaker attractions [iHSl- This is the background of the 
noted Weeks-Chandler- Andersen (WCA) papers [1, 0] in 
which a repulsive reference potential is arrived at by cut- 
ting off the potentials at their minima; similar notable 
works emphasizing the dominance of the repulsive inter- 
actions were published by Widom 0], Barker and Hen- 
derson 0, and Gubbins et al. During the 1970's 
it was demonstrated that the WCA reference represents 
well the liquid structure as monitored, e.g., by the radial 
distribution function 0, |3] • Since then thermodynamic 
perturbation theories [7[ start from the insight that in- 
termolecular attractions give rise merely to an almost 
constant mean-field background potential. 

There are several examples where the WCA refer- 
ence reproduces the liquid dynamics, but exceptions to 
this have also been reported Q. In a striking case, 
Berthier and Tarjus recently showed by simulation that 
the dynamics of the Kob-Andersen binary Lennard- Jones 
(KABLJ) viscous hquid is not at all properly reproduced 
by the WCA reference [£^]. Berthier and Tarjus concluded 
that, while the attractive forces have little effect on the 
liquid's structure, they affect the dynamics in a highly 
nontrivial and nonperturbative way. In view of these 
finding an obvious question is whether other repulsive 
potentials can better reproduce the KABLJ liquid's dy- 
namics. We focus here on inverse power-law (IPL) poten- 
tials, which are well known to have a number of simple 
properties We simulated three Kob-Andersen type 
binary mixtures: the original LJ version (KABLJ) [llj, 
the WCA version (KABWCA), and an inverse power-law 
version (KABIP). The KABLJ hquid is an 80:20 mixture 



of two LJ particles, A and B [11]. This liquid has become 
a standard system for studying viscous liquid dynamics, 
and like most binary mixtures with size ratio below 0.9 
it doesn't easily crystallize [13, [T^] ■ If Tij is the distance 
between particles i and j, the KABLJ potential energy 
is given by 



U 



(1) 



The indices i and j vary over all particles, and the en- 
ergy and length parameters are given as one of three 
combinations: eab/saa = 1-5 and Ebb/saa = 0.5, 
(tab/ctaa = 0.8 and (JbbIoaa = 0.88 [ll|. 

We wish to compare the dynamics of the KABLJ liq- 
uid to that of an IPL liquid. A simple IPL ansatz that 
respects the difference between A and B particles is 
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This ansatz leaves just two parameters to be determined. 



A and n. For an IPL liquid with 



interactions the 



virial obeys W — {n/3)U for aU microscopic states [lOj. 
Figure [Ija) shows constant-volume equilibrium fluctu- 
ations of virial versus potential energy of the KABLJ 
liquid at three state points with same density in the 
supercooled regime (the three colored ovals) [14|. The 
slope 7 of the best fit line is 5.16. Thus we choose 
n = 3 X 5.16 = 15.48 as the IPL exponent Once 
n is fixed, A is determined by plotting the fluctuations of 
the KABLJ potential energy versus J2i>j ^iji'^ij /''"ij)^ ^ 
referring to configurations drawn from the KABLJ sim- 
ulations (Fig. mb)). This plot yields A = 1.945 as the 
value ensuring overall agreement between the magnitudes 
of the KABLJ and the KABIP potential energy fiuctua- 
tions l™. 
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Kob-Andersen binary Lennard-Jones liquid ( KABLJ) 
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FIG. 1: (color online) (a) The IPL exponent is estimated 
to be n = 3 X 5.16 = 15.48 by evaluating the slope of a 
scatter plot of the virial and potential energy fluctuations of 
the KABLJ liquid at three state points with same density 
(open symbols, NVT simulations of 1000 particles). Average 
values of the state points are indicated by yellow x 's. (b) The 
IPL prefactor is estimated to be A = 1.945 from the slope 
in a scatter plot of ^^^^ ^ijif'ij / a-nd potential energy 
fluctuations. The number R gives the virial / potential energy 
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FIG. 2: (color online) Radial distribution functions of the 
AA, AB, and BB particle pairs of the KABLJ, KABIP, and 
KABWCA liquids at one state point. 



Before comparing the dynamics of the KABLJ and 
KABIP liquids we briefly consider their structure by eval- 
uating the radial distribution functions (RDF) in Fig. [5] 
This figure also includes the RDF of the KABWCA liq- 
uid. The IPL is slightly better than the WCA, but overall 
the three potentials give quite similar RDFs. The largest 
deviations are seen for the small particle (BB) distribu- 
tion; note that all three potentials capture its charac- 
teristic small first peak (with peak value below unity) 
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FIG. 3: (color online) Incoherent intermediate self-scattering 
function for the A particles of the KABLJ, KABIP, and 
KABWCA liquids along the p — 1.20 isochore for five tem- 
peratures. At the highest temperature the three models 
have the same dynamics, but at lower temperatures only 
the KABLJ and KABIP liquid have the same dynamics; the 
KABWCA liquid is considerably faster. Asterices (KAB- 
WCAIP) denote the IPL approximation to the KABWCA 
liquid (at each state point the parameters A and n are 
uniquely determined from the equilibrium virial and poten- 
tial energy fluctuations, resulting in at p = 1.20 for T — 
{0.45, 0.5, 0.6, 0.8, 2.0}: n = {20.67, 20.28, 19.59, 18.69, 16.38} 
and A = {1.176, 1.205, 1.231, 1.267, 1.471}). 



stemming from the strong AB affinity. 

Figure [3] compares the KABLJ, KABIP, and KAB- 
WCA dynamics via the A particle incoherent intermedi- 
ate self-scattering function evaluated at the wave vector 
of the first peak of the AA structure factor, q — 7.25 
(reduced units referring to the A particle parameters are 
used throughout). The simulations were performed over 
a range of temperatures on the p = 1.20 isochore, giving 
more than four decades of variation in the structural re- 
laxation time. The agreement between the KABLJ and 
KABIP dynamics is good. The intermediate scattering 
function of the KABWCA liquid likewise agrees with that 
of the KABLJ liquid at the highest temperature, but the 
KABWCA dynamics is too fast at low temperatures Q . 
Comparing Figs. [5] and [3] one concludes that the RDFs 
do not determine the dynamics. That pair correlation 
functions do not determine the dynamics may be rele- 
vant, for instance, for illuminating the applicability of 
mode-coupling theories [17j. As an aside. Fig. [3]includes 
also a system (KABWCAIP) approximating the KAB- 
WCA system by IPL potentials. As for the KABLJ hquid 
the two parameters A and n were determined (uniquely) 
from the equilibrium fluctuations (but in this case, inter- 
estingly, the parameters vary with state point whereas 
they are virtually constant for the KABLJ liquid). 

The structural relaxation time Tq is defined here as 
the time where the AA incoherent self-scattering func- 
tion is 1/e. This quantity's temperature dependence is 
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FIG. 4: (color online) (a) Structural relaxation time Ta as a 
function of inverse temperature defined as the time when the 
A A incoherent intermediate scattering function is 1/e at wave 
vector q = 7.25(p/1.2)^''^ . Lines between isochoric points are 
shown for clarity, (b) The reduced structural relaxation time 
of the KABLJ and KABIP liquids is a function of the variable 
p^/r, confirming density scaling 18, 19]. In both (a) and (b) 
the red symbols give the KABWCA results at densities 1.15, 
1.20, and 1.30. 



shown in Fig. HJ^a). We here varied the density between 
1.10 and 1.40, corresponding to pressures between -1 and 
11 kbar in the supercooled viscous regime (argon units). 
Again there is good agreement between the KABLJ and 
KABIP systems (Fig. SJa)). Figure ID^b) shows that the 
KABLJ liquid obeys density scaling [l8|,|l9|, i.e., that the 
reduced structural relaxation time is a function of p'^ /T . 
In particular, this implies that the IPL reference poten- 
tial also conforms well to the fragility scaling of Berthier 
and Tarjus (Fig. 2 in Ref. 11). 

The simulation results for Tq,(T) may be fitted by 
several well-known analytical functions for the temper- 
ature dependence along an isochore, both by functions 
with a dynamic divergence (e.g., the Vogel-Fulcher equa- 
tion) and by functions without a dynamic divergence 
(e.g., q uadratic scaling [20|] or other recent suggestions 
[2l|,|22|). For clarity of presentation we did not plot these 
functions. In all cases density scaling dictates the density 
dependence of the scaling parameters with the exponent 
7 determined from simulations at one density (Fig. [ij. 

Finally, we show in Fig. [5ja) the isochoric specific heat 
per particle, cy, as a function of temperature for the 
three potentials at different densities. Fig. [Sjb) shows 
that the KABLJ and the KABIP liquids both conform 
to density scaling as well as to Rosenfeld-Tarazona's cv 
scaling (cy-3/2 oc T-^/S) The KABLJ and KABIP 
liquids agree, while the KABWCA liquid gives somewhat 
smaller cy's. Thus not only for the dynamics, but also 




FIG. 5: (color online) (a) Isochoric specific heat per particle, 
Cv, in units of fcs. Lines between isochoric points are shown 
for clarity. The KABIP cy is closer to the KABLJ than is the 
KABWCA. (b) The same data scaled according to the density 
scaling recipe in conjunction with the Rosenfeld-Tarazona cv 
scaling [2^ , which is indicated analytically by the blue dashed 
line (the number 1.25 being a fitting parameter). 



for a static thermodynamic quantity like cy does an IPL 
system provide a superior reference system. 

Why does an IPL potential so well reproduce the 
KABLJ liquid's dynamics? It was recently argued that 
van der Waals-type liquids have a "hidden" (approxi- 
mate) scale invariance, which is reflected in the fact that 
these liquids to a good approximation inherit a number 
of properties of repulsive inverse power law (IPL) poten- 
tials As mentioned, an IPL potential exhibits 100% 
correlation between virial W and potential energy U for 
all microscopic configurations \W = (n/3)[/]. Hidden 
scale invariance is manifested in the appearance of strong 
correlations of the thermal equilibrium WU fluctuations 
in the NVT ensemble (Fig. [iJ; it comes about as fol- 
lows 



24|. Highly asymmetric potentials like the LJ pair 



potential are well approximated around the minimum 
by an "extended inverse power law" (eIPL) of the form 
v{r) — Ar~" + B + Cr. When summing over all nearest- 
neighbor distances, the linear term almost doesn't fluc- 
tuate if the total volume is fixed. Consequently, several 
properties of strongly correlating liquids are well repro- 
duced by IPL potentials. It should be noted, though, 
that because the linear term of the eIPL potential gives 
a volume-dependent contribution to the free energy, the 
IPL free energy, pressure, compressibility, and equation 
of state, represent poorly the liquid [2^ . 

Strongly correlating liquids include the van der 
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Waals and metallic liquids, but neither covalently and 
hydrogen-bonded liquids, nor strongly ionic liquids 24 1 
- this is because competing interactions generally spoil 
strong WU correlations. To a good approximation 
strongly correlating liquids obey density scaling [18] ac- 
cording to which the relaxation time at varying tem- 
peratures and pressures is a function of p'^ /T. In sim- 
ulations of strongly correlating liquids 7 is determined 
from the approximate identity AW{t) = jAU{t) (Fig. 
H]) O 24 1 . Generally, strongly correlating liquids are 
characterized by the existence of "isomorphic" curves in 
their phase diagram along which several quantities are 
almost invariant, e.g., the excess entropy and the relax- 
ation time in reduced units (25l |. Whenever 7 is fairly 
constant throughout large parts of the phase diagram - 
as appears to be the case for most strongly correlating 
liquids [IJl - an isomorph is to a good approximation 
given by p"' /T = Const. This explains density scaling for 
the class of strongly correlating liquids [25 1. 

In conclusion, structure, dynamics, and some thermo- 
dynamic properties of the KABLJ liquid are well repre- 
sented by a repulsive IPL reference. Our findings suggest 
that the WCA reference fails to reproduce the viscous dy- 
namics because this reference does not incorporate the 
hidden scale invariance of the KABLJ liquid. 
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